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ABSTRACT 

Hard X-ray surveys have proven remarkably efficient in detecting intermediate po- 
lars and asynchronous polars, two of the rarest type of cataclysmic variable (CV). 
Here we present a global study of hard X-ray selected intermediate polars and asyn- 
chronous polars, focusing particularly on the link between hard X-ray properties and 
spin/orbital periods. To this end, we first construct a new sample of these objects by 
cross-correlating candidate sources detected in INTEGRAL/IBIS observations against 
catalogues of known CVs. We find 23 cataclysmic variable matches, and also present 
an additional 9 (of which 3 are definite) likely magnetic cataclysmic variables (mCVs) 
identified by others through optical follow-ups of IBIS detections. We also include in 
our analysis hard X-ray observations from Swift /BAT and SUZAKU /HXD in order 
to make our study more complete. We find that most hard X-ray detected mCVs have 
Pspin 1 1 Porb < 0.1 above the period gap. In this respect we also point out the very low 
number of detected systems in any band between P 8p in/ Porb = 0.3 and P sp in/Porb = 1 
and the apparent peak of the P sp i n /Porb distribution at about 0.1. The observational 
features of the P sp i n — Porb plane are discussed in the context of mCV evolution 
scenarios. We also present for the first time evidence for correlations between hard 
X-ray spectral hardness and P sp im Porb and P sp in/ Porb- An attempt to explain the 
observed correlations is made in the context of mCV evolution and accretion footprint 
geometries on the white dwarf surface. 

Key words: gamma-rays: observations, surveys, X-ray binaries, magnetic cataclysmic 
variables 



1 INTRODUCTION 

Cataclysmic variables (CVs) are close binary systems con- 
sisting of a late-type star transferring material onto a white 
dwarf (WD) via Roche-lobe overflow. Magnetic CVs (mCVs) 
are a small subset of the catalogued CV s (« 10% — 20%, 
iDownes et al. I l2005l : lRitter&~ Kolb 2003), and fall into two 
categories: polars (or AM Her types after the prototype sys- 
tem) and intermediate polars (IPs or DQ Her types). The 
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WDs in polars possess such strong magnetic fields that they 
can synchronise the whole system, yielding P or t = Pspin- 
The strong magnetic field in these systems is confirmed by 
strong optical polari sation and measurments of cyclotron 
humps (|Warnerll200&t ). Accretion in polars is thought to fol- 
low the magnetic field lines of the WD straight from the LI 
point onto the WD magnetic poles, and no accret ion disk is 
expected (for a review of polars, see ICropperlll99(ih . Another 
possible class of systems called asynchronous polars (APs) 
might exist, where the spin and orbital periods are out of 
synchronisation by only a few percent. It is not known ex- 
actly why this is but one suggestion is that these systems 
are polars which have had a rec ent nova event , kicking them 
slightly out of synchronisation l|Warneri r2003h For IPs, the 
lack of strong optical polarisation implies a weaker magnetic 
field, not powerful enough to synch ronise the secondary (for 
a review of IPs, see |Patterso4l993h . In these systems, mate- 
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rial leaving the LI point usually forms an accretion disc up 
to the point where the magnetic pressure exceeds the ram 
pressure of the accreting gas. From this point onwards the 
accretion dynamics are governed by the magnetic field lines, 
which channel the material onto the WD magnetic poles. 
The nature of these systems is confirmed by the detection 
of coherent X-ray modulations associated with the spin pe- 
riod of the WD. 

In the simplest scenario for X-ray production in 
mCVs, the magnetically channeled accretion column im- 
pacts the WD poles producing hard X-rays from thermal 
bremsstrahlung cooling by free electrons with kT of the or- 
der of 10s of keV !|Warnerl2003l ; |Cropperlll990h ■ The emission 
is thought to originate in the post-shock region, a region be- 
low the shock front created from the impacting accretion 
column. Softer X-rays are also produced from the absorp- 
tion and reprocessing of these higher energy photons in the 
WD photosphere. As a result, both polars and IPs are ex- 
pected to emit high energy photons, but discrepancies exist 
between the observed ratio of soft-to-hard X-rays between 
polars and IPs, with p o lars showing an excess of observed 
soft X-rays (|Lamblll98a ). IChanmugam et al.l (| 1991T ) and oth- 
ers have reported that the total X-ray luminosities of IPs are 
greater than those of polars by a factor of ~ 10, attributed 
mainly to the higher accretion rates. Moreover it has been 
proposed that strong magnetic fields in pola rs produce a 
more "blobby" flow than in IPs (|Warnerll2003r ). These high 
density "blobs" are then able to penetrate within the post- 
shock region, emitting fewer bremsstrahlung photons and 
contributing more to the observed X-ray blackbody spec- 
tral component, thought to be produced at the base of the 
post-shock region. 

In recent years an increasing number of mCVs have 
been detected and dis covered by hard X-ray telescopes suc h 
as INTEGR AL/IBIS dBarlow et al.l2006l ; ILandi et aill2009l ). 
Swift/BAT (IBrunschweiger et alJl2o69T Tand SUZAKU/HXD 
<|Teradaet al.ll200Sl ). Here, we present a newly compiled 
list of IBIS mCVs, forming part of a larger compilation of 
hard X-ray selected IPs from all three telescopes mentioned 
above, allowing us to study the general properties of this 
hard X-ray selected sample. 

This selected sample of CVs will allow us to investi- 
gate the observational properties of the P or b~Pspin plane, 
and compare them to mCV evolution models. Moreover, we 
will investigate the spectral properties of our sample, and 
present for the first time observed spectral hardness correla- 
tions observed within IPs as a function of their orbital and 
spin parameters. 



2 RECENT HARD X-RAY OBSERVATIONS 

The INTEGRAL satellite, launched in October 2002, has 
now carried out more than 6 years of observations in the 
energy range 5 keV to 10 MeV. In particular, the INTE- 
GRAL/IBIS surve y is one of the main mission objectives . 
The IBIS detector (|Ubertini et alj|2003t iLebrun et al.ll2003l ) 
has been optimised for survey work, with a large field of 
view (30°) and with unprecedented sensitivity in the soft- 
gamma ray regime, yielding excellent imaging capabilities. 
It is worth emphasising that the IBIS survey has been op- 
timised to detect faint persistent sources, which mCVs are. 



The aim of the survey is to expand the current knowledge of 
the 20-100 keV sky by cataloguing high-energy sources and 
examining their properties, both individually and globally. 
The IBIS survey dataset consists of dedicated observations 
along the Galactic plane and around the Galactic centre. 
Additionally, a combination of pointed and deep exposure 
observations are added to the dataset once they become pub- 
lic. As a result, the latest release of the INTEGRAL /IBIS 
survey provides all-sky coverage, albeit with spatially vari- 
able sensitivity. The depth of t he IBIS survey has increased 
significantly with each release (|Bird et alj|2004 l2006t l2007f ) 
and has now reached a peak sensitivity corresponding to a 
flux limit below 1 mCrab in th e 20-100 keV range . The latest 
study on INTEGRAL CVs by ILandi etlll fcOOSft contained 
22 sources found within the IBIS /ISGRI survey or identified 
as CVs later on through optical foll ow-up. 

This work com plements that of ILandi et alj (120091 ) and 
iBarlow etT al, (2006), in particular including a larger dataset 
yielding deeper observations. The IBIS dataset used con- 
sists of over 36,000 individual Science Window (ScW) point- 
ings, covering 6 years of observations (revs w 46-660). These 
have been processed with the latest pipeline ( OS A 7. 0, 
iGoldwurm et al.ll2003l ). and mosaics were created from the 
deconvolved images in 5 energy ranges 1 . Staring and per- 
formance verification (PV) observations have not been used 
for the mosaic creation, and noisy ScWs have been excluded 
based on the RMS of the individual images. These new sur- 
vey maps are a considerable improvement on any previous 
ones, due to the new pipeline software, together with the 
significantly increased exposure times. 

Swift, has been optimised to locate gamma-ray bursts, 
and as a consequence the main hard X-ray instrument, the 
Burst Alert Monitor (BAT), has similar capabilities to IBIS, 
possessing a large field of view and operating in essentially 
the same energy range. BAT has also been used for sur- 
vey work ( Tu eller et al.l 120091 ) , and in particular has als o 
detected a high number of IPs (|Brunschweiger et al.l [20091 ) . 
In order to make our study as c omplete as possible, we have 
decided to include the IPs from IBrunschweiger et al] (120091 ) 
observed with the Swift/BAT detector. Moreover, one extra 
IP has been added, AE Aqr, as observed by SUZAKU/HXD 
ijTerada et al.ll200a ). sampling again a very similar energy 
range to IBIS. The similarity of energy range allows us to 
construct a hard X-ray selected sample with minimal biases. 

Prior to this work, hard X-ray detections have been re- 
ported for ss 30 mCVs above 10 keV. More than 90% of 
these are IPs, and there are also two rare asynchronous po- 
lars. When compared to the older soft X-ray selected sam- 
ples of mCVs the picture is somewhat different. First of all, 
as one would expect, soft X-ray detectors are more sensi- 
tive to mCVs and consequently will produce larger samples; 
including equal amounts of polars and IPs. 



3 THE HARD X-RAY CV POPULATION 
3.1 INTEGRAL/IBIS CVs 

Here we describe the catalogue matching procedure adopted 
in order to identify known CVs in our IBIS dataset. To do 

1 20-40keV, 30-60keV, 20-100keV, 17-30keV and 18-60keV 
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this we require the most complete and up-to-date catalogue 
of such objects, thus we merge the two most complete CV 
catalogues in the lit erature: The Catalog ue and Atlas of Cat- 
aclysmic Variables (|Downes et al.ll2005l . here after DWScat) 
and t he Catalogue of Cataclysmic Binaries l|Ritter fc Kolbl 
12001 hereafter RKcat). DWScat contains 1830 CVs whilst 
RKcat contains 731, and we note that 656 CVs are common 
to both catalogues. The principal reason for RKcat having 
fewer objects is that only CVs with known orbital periods 
are included in the sample; however, RKcat also includes a 
few CVs that DWScat does not report. Our final known CV 
set therefore contains 1905 CVs (hereafter DRKcat). Fewer 
than w 9% of the total number of CVs within DRKcat are 
known to be magnetic in nature (included in the catalogue 
as either DQ Her, AM Her or IP), and only approximately 
3% (56 sources) are IPs. 

Catalogue matching has been performed between the 
total CV set produced, which contains 1905 CVs and a pre- 
liminary IBIS candidate excess list containing real sources 
constructed in the same way as for the ISINA algorithm 
(|Scaringi et al ]|2008l . Section 3), but now on a larger dataset, 
containing over 9000 excesses constructed from the public 
data available for revolutions 46 -660. This includes excesses 
detected in the final mosaics, revolution mosaics and revo- 
lution sequence mosaics in any of the 5 main energy bands 
in order to locate variable candidates too. 

W e have carrie d out the same correlation exercise 
as in Barlow et al. (2006), usin g a similar method to 



IStephen et "all (|2006l '). Similarly to iBarlow et al.1 l|2006h we 



employ a search radius of 4', which also corres ponds well 
with the expected error on faint IBIS detections (|Gros et al.l 
2003). For a search radius of 4' we obtain 56 sources as con- 
firmed or candidate CVs, of which 23 are expected to be 
false coincidences. We have visually inspected all of the cor- 
related sources and found 33 matches coincide with mainly 
non-CV globular cluster sources and previously identified X- 
ray objects; however some are image artefacts related to the 
Galactic centre region. It is important to note that with a 
~ 2' typical source location accuracy it is very hard to asso- 
ciate a detection with an optical counterpart alone. We have 
performed the same exercise by increasing the search radius 
to 5' which increases the sample to inspect to 76 candidates. 
We find that all are false matches with the possible excep- 
tion of the Dwarf Nova DN V1830 Sgr located 4.8' away 
from the IBIS detection detected in revolution 106 (MJD 
53128.8 - 53131.7). This is slightly outside the 90% error ra- 
dius for a ~ 6<t detection and we cannot definitely associate 
the two at the moment. 

Table Q] shows the main characteristics of the 23 ob- 
jects identified from our correlation analysis. We have es- 
timated di stances f or th is sample using the method de- 
scribed by iKniggel <|2006h . based on the evolution, of the 
dono r star where 2MA SS K-band magnitudes were avail- 
able (|Cutri et al.ll2003h and Porb < 6.2 hours 2 . In addition 
we show in Table [2] the 9 other IBIS-detected mCVs used in 
this work. These were not part of the correlation analysis, 
because they were not present in DRKcat, but have been 
identified through optical spectroscopy following the IBIS 
discovery. Of the 23 objects considered to be real matches 
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from our analysis, 17 are previously known INTEGRAL de- 
tected CVs, whilst 6 sources are new detections. Most of 
the new objects are of the intermediate polar subclass with 
the possible exception of TW Pic, which is cons idered by 
some as a possible VY Scl star (|Norton et al.ll20o"oT ) . and AX 
J1832. 3—0840, which is not identified in full at the moment. 

In total this paper presents 32 IBIS-detected CVs, 
increasing the nu mber from the latest previous study of 
lLandi et all (|2009T ) by 10 systems. 



3.2 Swift/BAT and SUZAKU/HXI) CVs 

Swift /B AT has also observed a large number of mCVs, 
and it is worth including these in our study for complete- 
ness. We decided to inclu de all the 22 BAT detected IPs 
l|Brunschweiger et al.ll20o"9T ). where 14 have been observed by 
IBIS as well. Similarly to the IBIS-only IPs, the BAT-only 
IPs are all placed above the period gap with the exception 
of EX Hya. 

SUZAKU /HXD on the other hand has a different ob- 
serving strategy composed of small field of view pointings. 
This does not allow the telescope to produce survey data 
like IBIS or BAT, however mCVs have been detected and 
observed wit h HXD. In particul ar HXD had observed the 
IP AE Aqr l|Terada et al.l l2008h . which has not been ob- 
served with either IBIS or BAT, and therefore is included 
in our study as well. We caution however that the possible 
origin of the hard X -rays in AE Aqr could be non-thermal 
l|Terada et al.l 12008) . This however has not been shown in 
full, and given the power-law model fit to the 3-25 keV 
spect ra of AE Aqr yielding an index of 2.10 (|Terada et al.l 
2008) not far from the indeces found in fast spinning mCVs 
l Landi et al.ll2009h we believe this object should still be in- 
cluded in our analysis. 

All the IPs observed with BAT and HXD used in this 
study are presented in Table [3] together with their orbital 
and spin periods. 



4 HARD X-RAY SELECTED MAGNETIC CV 
PROPERTIES 



As revealed in lBarlow et all (120061 ) and lLandi et al.l (2009), 
the vast majority of IBIS detected CVs are magnetic (the 
only exception being SS Cyg). The total number of objects 
with known spin and orbital periods in our sample is 30 
systems, of which 22 are IBIS detections (18 IPs). It is in- 
teresting to note the relative incidence of polar systems and 
intermediate polars. Optically selected samples favour the 
former, with the number of known polars being twice that 
of IPs. However, in our hard X-ray selected sample, the pic- 
ture is very different, with only 4 polars being included in 
the sample (2 of which are APs). This result is expected, 
as IPs are known to produce « 10 times more hard X-rays 
than polars due to the ir higher mas s transfer and intrinsi- 
cally harder spectrum |Warncr 2003|). 

Only two synchronous polars have been detected in our 
sample, both at relatively close distance and in regions of the 
sky with high exposure times. This leads us to conclude that 
the flux in the hard X-ray range for these systems is much 
lower when compared to IPs or APs. We note from Table [1] 
that the two detected polars are among the closest objects 



Table 1: Results of the IBIS — DRK catalogue matching with 4' search. 
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a,8 b 
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Porb 


p 
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Refs 




(IBIS position) 
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(ct s" 1 ) 


(ks) 


20-100 keV 


(min) 


(s) 


(PC) 




1RXS J002258. 3+614111 


5 739 61 71 4 


IP 


1.7 


B4(8.7) 


15 + 001 


3767 


0.81 


241.98 


563 53 


510 


n 2 3 41 


V709 Cas 


7.207,59.303 


IP 


0.8 


B5(54.3) 


1.03 ± 0.01 


3562 


5.53 


320.4 


312.77 


300 


\1 2 51 


XY Ari* 


44 047 1 9 457 


IP 


1.1 


R5C5 B) 


53 + 12 


119 


2.85 


363 884 

*-j \j '.J . t_j rt 


206 298 


610 


M 2 61 


GK Per 


52 777 43 928 


IP /DN 


1.7 


B5(4.7) 


26 + 07 


277 


1.4 


2875.4 


351.34 




[1 21 


TV Col* 


82 357 -32 81 9 


IP 


0.1 


B4(11.6) 


68 + 08 


248 


0.37 


329.181 


1911 


330 


[1,2,7,8] 


TW Pic* 


83 766 -57 998 


IP?/VY Scl?»» 


2.5 


Bl(5.8) 


3 + 07 

yj ' *-J _1_ yj > \J t 


363 


1.61 








[i 2 9 10 111 


BY Cam 


85 728 60 842 


AP 


1.2 


B5(5.1) 


69 + 11 

yj . \j *j _l yj . -i- J. 


162 


3.7 


201.298 


11846.4 


140 


\1 21 


MU Cam 


96.316,73.567 


IP 


0.6 


B4(5.4) 


0.24 ± 0.06 


548 


1.29 


283.104 


1187.24 


440 


[1,2,12] 


SWIFT 10732 5— 1331* 


113.13 -13.513 


IP 


1.6 


B3(6.1) 


39 + 06 


409 


2.09 


336.24 


512.42 




[1,2,13] 


V834 Cen 


212.260,-45.290 


P 


0.9 


Bl(5.4) 


0.16 ± 0.03 


1675 


0.86 


101.51712 


6091.0272 


70 


[1,2] 


ICR J14536-5522 


223.421,-55.394 


P 


2.0 


B4 (11.9) 


0.27 + 0.03 


2658 


1.45 


189.36 


11361.6 


110 


[1,2,14] 


NY Lup 


237.052,-45.481 


IP 


0.5 


B5(49.1) 


1.17 + 0.03 


3141 


6.28 


591.84 


693.01 




[1,2,15] 


V2400 Oph 


258.173,-24.279 


IP 


2.2 


Bl(33.4) 


0.68 ± 0.02 


4453 


3.65 


204.48 


927.6 


180 


[1,2,16] 


1H 1726-058 


262.606,-5.984 


IP 


0.7 


B5(22.8) 


0.85 ± 0.04 


1449 


4.56 


925.27 


128 




[1,2,17] 


V2487 Oph 


262.960,-19.244 


IP/Nm 


2.3 


B3(9.1) 


0.18 ± 0.02 


4562 


0.97 








[1,2,18] 


AX J1832.3-0840* 


278.083,-8.721 


? 


3.1 


B4(5.5) 


0.07 + 0.03 


3090 


0.38 








[1,2,19,20] 


V1223 Sgr 


283.753,-31.153 


IP 


0.8 


B5(52.2) 


1.45 ± 0.03 


2358 


7.79 


201.951 


746 


150 


[1,2] 


V1432 Aql 


295.052,-10.421 


AP 


0.2 


B5(10.8) 


0.69 ± 0.07 


429 


3.7 


201.938 


12150.4 


240 


[1,2,21,22] 


V2069 Cyg 


320.906,42.279 


IP. 


1.8 


B5(6.2) 


0.21 ± 0.03 


1648 


1.13 


448.824 


743.2 




[1,2,23,24] 


1RXS J213344. 1+510725 


323.446,51.122 


IP 


0.3 


B5(25.8) 


0.65 ± 0.03 


2207 


3.49 


431.568 


570.82 




[1,2,25] 


SS Cyg 


325.698,43.582 


DN 


0.6 


B5(23.0) 


0.7 + 0.03 


1674 


3.76 


396.1872 






[1,2,26] 


FO Aqr 


334.514,-8.354 


IP 


1.7 


B4(6.1) 


0.65 + 0.2 


54 


3.49 


290.966 


1254.45 


250 


[1,2,27] 


AO Psc* 


343.815,-3.194 


11' 


1.3 


B4(4.8) 


0.43 ± 0.11 


108 


2.31 


215.461 


805.2 


200 


[1,2,28,29] 



[7 
[11 
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. a * indicates new hard X-ray detections 
^ Right ascension and declination in degrees, J2000 

~!S c IP=intermediate polar, P=polar, AP=asynchronous polar, N=nova, DN=dwarf nova. All are confirmed except for »:probable, ••: possible 

'g^ ^ Angular distance between the DRKcat catalogue positions and the IBIS coordinate 

S e IBIS detection only. Map with maximum significance: (Bl) 20-40 keV; (B2) 30-60 keV; (B3) 20-100 keV; (B4) 17-30 keV; (B5) 18-60 keV; in brackets appears the significance value. 

^ J Determined between 20-100 keV 

9 The flux is calculated assuming a power law spectra with index of 2.9, the average index for IPs fearlow et al,|[2006l) 

^ The distances have been computed with 2MASS K band magnitudes (Cutri ct al. 2003J using the method presented by[Kniggc (20Q6|) based on the evolution of the secondary. 

to 



Table 2: Additional mCVs detected by IBIS not included in DRKcat. 



Name 


a,8 a 


type" 


Map code c 


Count rate * 


Exposure 


Flux e 


Porb 


p 

1 spin 


Distanced 


Refs 










(ct s- 1 ) 


(ks) 


20-100 kcV 


(min) 


(k) 


(pc) 




IGR J08390-4833 


129.705,-48.524 


IP». 


B5(6.3) 


0.08 ±0.02 


3072 


0.43 








[1] 


XSS J12270-48599 


187.004,-48.906 


IP.. 


B5(9.9) 


0.42 ± 0.04 


955 


2.26 








[2,3,4] 


IGR J15094-6649 


227.406,-66.823 


IP 


B5(11.0) 


0.18 ±0.03 


3265 


0.97 


353.40 


809.42 


600 


[3,5] 


IGR J16167-4957 


244.140,-49.974 


IP.. 


B4(20.7) 


0.4 ±0.02 


3466 


2.15 








[3] 


IGR J16500-3307 


252.491,-33.114 


IP/DN 


Bl(13.3) 


0.33 ±0.03 


3353 


1.77 


217.02 


597.92 


270 


[5,6] 


IGR J17195-4100 


259.906,-40.997 


IP 


B5(23.3) 


0.54 ±0.02 


4279 


2.9 


240.30 


1139.55 


120 


[2,3,5] 


IGR J18173-2509 


274.353,-25.158 


IP.. 


Bl(14.9) 


0.27 ±0.01 


5744 


1.45 








[7] 


IGR J18308-1232 


277.696,-12.532 


IP.. 


B5(7.1) 


0.18 ±0.03 


3265 


0.97 








[8] 


IGR J19267+1325 


291.670,13.425 


IP. 


B5(7.5) 


0.15 ±0.03 


2963 


0.81 








[9,10] 



A ll classifications ar e co rrect except for , .: probable classificat ion ..: possible c lassification. References: [1 Kniazcv ct al. (2008); [2 Butters ct al. (2008) 
[4 lSaitou etall d2009l ); [5 lPretorius1 i l2009h ; [6 ]Masetti et al.l [120081 ); f7 ]Masetti et al.l d2008l) ; [8 lParisi et al] i2008h ; [9 lSteeghs et al.l J2OO8I I ; [lO lEvans et alj J20Q. 



a Right ascension and declination in degrees, J2000 

^ IP=intermediate polar, P=polar, AP=asynchronous polar, N=nova, DN=dwarf nova. All are confirmed except for «:probable, ..: possible 

c IBIS detection only. Map with maximum significance: (Bl) 20-40 keV; (B2) 30-60 keV; (B3) 20-100 keV; (B4) 17-30 keV; (B5) 18-60 keV; in brackets appears the significance value. 
d Determined between 20-100 keV 

e The flux is calculated assuming a power law spectra with index of 2.9, the average index for IPs l lBarlow et al.ll2006l ) 

/ The distances have been computed with 2MASS K band magnitudes dCutri et al.ll2003li using the method presented by |Kniggc ( 200(|) based on the evolution of the secondary. 
9 Note that there are conflicting results regarding the nature of this object. RXTE data suggests an IP nature, whilst SUZAKU data suggests an LMXB nature. See references. 
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Table 3. Additional IPs used in this work detected with Swift /BAT and SUZAKU /HXD with known spin and orbital periods. 



Name 


a, S 


type 


Detection 


Porb 


p 

1 spin 


Refs 










(min) 


w 




V1062 Tau 


75.615,24.756 


IP 


BAT 


597.60 


3726 


[1] 


TX Col 


85.834,-41.032 


IP 


BAT 


343.15 


1909.7 


[1] 


V405 Aur 


89.897,53.896 


IP 


BAT 


249.12 


545.455 


[1] 


BG CMi 


112.871,9.940 


IP 


BAT 


194.04 


847.03 


[1] 


PQ Gem 


117.822,14.740 


IP 


BAT 


310.80 


833.4 


[1] 


DO Dra 


175.910,71.689 


IP 


BAT 


238.139 


529.31 


[1] 


EX Hya 


193.102,-29.249 


IP 


BAT 


98.257 


4021.62 


[1] 


AE Aqr 


310.038,-0.871 


IP 


HXD 


592.785 


33.0767 


[2] 



References: [l'Brunschwcigcr et al. (2009); [2'Terada ct al. (2008). 



in our sample, a nd hence it is prob ably for this reason that 
IBIS (and BAT. iTueller et aljlioosh can see them. We have 
also checked this by inspecting where other polars sit in the 
INTEGRAL exposure map and conclude that more deep 
observations are required before more of these systems are 
detected in the hard X-ray range. 

Two of the four confirmed asynchronous polars (APs) 
are also included in our sample. At first one might think 
that these systems should have properties resembling t he po- 
lar class. However, as shown bv lSchwarz et alj (|2005f ) . their 
accretion rates are ~ 10 — 20 times greater than that of 
polars. Moreover both INTEGRAL detected APs are 2-3 
magnitudes brighter than the two non-detected APs in the 
K-band. This, together with our distance estimates suggests 
that our efficiency at detecting APs in the hard X-ray range 
is w 50%, similar to the IPs. 

4.1 The P rb - Pspin plane 

Theoretical si mulations on the evol u tion o f mCVs have been 
performed by iNorton et al . (200§j, 120041 ). One interesting 
prediction of these models is that of different accretion flows 
for IPs depending on where they are in their evolutionary 
stage. In particular these models also predict the existence 
of spin equilibria among IPs, one at about P ap in/Porb ~ 0.1 
and a second one at P ap in/Porb ~ 0.6 (depending on mass 
ratio) . 

Figure[T]shows the Porb — Pspin plane for mCVs together 
with their orbital, spin and synchronicity (Pspin/Porb) distri- 
butions. The hard X-ray detected systems used in this work 
are shown with stars. Also plotted for reference is the period 
gap and three "synchronicity" lines showing Pspin = Porb 
(polars), Pspin = 0.1P or 6 and P sp in = 0.3P or6 . It is clear 
that most of the hard X-ray detected systems are above the 
period gap and have P sp i n ^ 0.1P or ;,. The only IP system 
outside this range is EX Hy a which is closer than 100 parsec 
(jBrunschweiger et al.ll2009t ). Similarly to the two detected 
polars, it is of no real surprise that hard X-ray detectors 
can see this close object. 

The fact that no IP has yet been observed with hard 
X-ray telescopes above the period gap with P sp in/ Porb > 0.1 
suggests that these IPs have different accretion flows com- 
pared to IPs with Pspin / Porb < 0.1. This idea is supported 
by the fact that the distribution of all known mCVs does in- 
deed seem to peak at about P sp in/Porb ~ 0.1 in the bottom 
left panel in Fig[TJ where a spin equilibrium has been pre- 
dicted. As IPs evolve from low synchronicity (Pspin/Porb « 
0.1) towards their Pspin/Porb ~ 0.1 equilibrium, their accre- 



tion flows resemble those of propeller systems, where a lot 
of the material incoming from the secondary is actually pro- 
pelled away and does not reach the pole of the WD. So in 
essence most IPs with Pspin / Porb « 0.1 have yet to reach 
their equilibrium (regardless of field strength), and their ac- 
cretion flows are differ ent from any other IP e lsewhere in 
the Pspin - Porb plane j|Norton et al.ll2008l . [200i ). As a con- 
sequence one would not necessarily expect the hard X-ray 
emission mechanisms to be the same for all IPs. 

Another interesting observational feature of this plane 
is that only one unconfirmed IP system, V697 Sco, lies 
within what we call the synchronicity gap: a region in the 
Porb — Pspin plane within Pspin > 0.3P or b and total syn- 
chronicity above the period gap. The low number of IP 
systems with high synchronicity above the period gap is 
also p r edicte d by the models. As explained bv lNorton et al.1 
(2008, 2004), as mCVs evolve, both their mass ratios and 
orbital period decrease. These trends individually cause op- 
posite shifts in the spin-to-orbital ratio at which the spin- 
equilibrium occurs for a given magnetic field. As a result, 
typical IPs with magnetic field strength of a few MG will 
evolve from being disc-like accretors at long orbital period 
(where Pspin / Porb ~ 0.1), to ring-like accretors at short or- 
bital period (where Pspin/Porb ~ 0.6), providing that they 
do not synchronise along the way and become polars. If the 
magnetic field of the systems is of the order of a few hun- 
dred MG, on the other hand, then they will be prone to 
synchronise and become polars, crossing the synchronicity 
gap fairly quickly, helping explain the low number of sys- 
tems in this region. This is the likely history of EX Hya and 
systems neighboring in the P orb — Pspin plane. 

We therefore predict that, given the already long expo- 
sure times accumulated with IBIS and BAT, the IPs within 
Pspin = O.lPorb and Pspin = 0.3P or b region will not be de- 
tected in significant numbers. Conversely we expect most of 
the IPs with Pspin < 0.1P or b to be observable with longer 
exposure times and better sensitivity. We also do not expect 
many systems to exist within the synchronicity gap. 

4.2 Are hard X-ray detected mCVs different? 

It should be clear by now that hard X-ray telescopes are 
more sensitive at detecting the IP population rather than the 
polar one. However it is not yet clear whether hard X-ray 
telescopes are producing populations which are consistent 
with being drawn from the general mCV population. In or- 
der to assess this we have performed a Kolmogorov-Smirnov 
test (KS test) on all the Pspin, Porb and Pspin/Porb distribu- 



Hard X-ray mCVs 7 




Figure 1. P or b vs. P sp in for mCVs taken from RKcat. Stars indicate mCVs detected at hard X-ray energies. Also plotted is the period 
gap and "synchronicity" lines. For reference we also display the orbital, spin and synchronicity distributions. The shaded green areas 
represent hard X-ray selected mCVs. 



tions of hard X-ray selected samples versus the known mCV 
population taken from RKcat. In all cases the test rejects 
the null hypothesis that the distributions are drawn from 
the same parent with > 99.99% confidence. 

As mentioned before, the difference in distributions 
within Port, (bottom panel in Fig.HJ between these sets could 
be caused by the fact that all mCVs below the period gap 
are intrinsically less luminous given their lower accretion 
rates. However it does not exclude the possibility that the 
X-ray emission mechanisms for mCVs below the period gap 
is substantially different than the mCVs above the gap. Both 
of these effects can be regarded as systematic, however the 
possibility of hard X-ray missions detectin g homogeneous 
IP sam ples has already been suggested by iGansicke et al.l 
(2005), and here we bring forward more observational evi- 
dence for this. It is clear however that if mCVs below the 
period gap emit hard X-rays, then current telescopes are not 
sensitive enough to detect them. 

Moreover, from the P sp in/Porb KS test (bottom left 
panel in Fig. [1} result we can also comment on the fact that 
hard X-ray telescopes are not sensitive to high synchronic- 
ity systems. These do not necessarily have to live below the 
period gap, and in fact about half of the polar population 
lives above the gap. This suggests that hard X-ray surveys 
are very insensitive to the polar population as well as mCVs 
below the period gap. We can also look at this result and 



suggest that, on the other hand, hard X-ray surveys are ex- 
tremely sensitive to IPs with long P or b (i.e. above the period 
gap) and to systems with low P S pin/ Porb- 

It is hard to say whether mCVs below the period gap 
and polars in general have different emission mechanisms 
with a KS test alone. However we believe that the evidence 
arising from the orbital and synchronicity distributions sug- 
gest this very strongly. 



5 HARDNESS PLANE CORRELATIONS 

The production of hard X-ray photons from mCVs is 
thought to originate in the post-shock region of the WD 
by bremsstrahlung cooling of free electrons. This is some- 
what different to softer X-rays (< 2 keV) seen from mCVs, 
which can originate from a blackbody component close to 
the WD surface. In fact in recent years many medium 
resolution X-ray spectra have been o btained for di f ferent 



Schwarz et al.l 20051; iButters et all 120081; 


Evans & Hellierl 


2007: bone & Maedziarz 1998: Landi et al. 


2009) and have 



been fitted with a soft blackbody component (kT « 80eV) 
plus a hard comp onent character i sed b y the stratified ac- 
cretion column of ICropper et all (| 19991 ). For those mCVs 
that are detected in the hard X-ray range the detection only 
samples the bremsstrahlung component. In particular the 
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Figure 2. 30-60/17-30 koV hardness versus spin period for the 
hard X-ray selected mCVs used in this work. Polars and APs are 
shown in empty squares. 
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Figure 4. 30-60/17-30 keV hardness versus Pspin/Porb for the 
hard X-ray selected mCVs used in this work. Polars and APs are 
shown in empty squares. We note the evident correlation for IPs 
With Pspin/Porb < 0.1 
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Figure 3. 30-60/17-30 keV hardness versus orbital period for the 
hard X-ray selected mCVs used in this work. Polars and APs are 
shown in empty squares. 



hard X-ray energy range (> 17 keV) is telling us about the 
temperature distribution of components within the multi- 
temperature bremsstrahlung emission, not the ratio of hard- 
to-soft X-ray components. 

Keeping this in mind, Fig. [2] [3] and [4] shows the scatter 
plots for hardness, defined as the count ratio in the 30-60 

keV and 17-30 keV bands, VerSUS Pspin , Porb and Pspin /Porb 

respectively for all hard X-ray detected mCVs used in this 
work. In red, we show all mCVs seen by IBIS, in blue, BAT- 
detected mCVs and in black, we show the only IP used 
in this work observed by HXD, AE Aqr. In order to ob- 
tain hardness ratios for BAT and HXD mCVs, we have re- 
produced their bremsstrahlung spectra (power law for AE 
Aqr) using the temperatures (or photon index ) prov ided 



A qr) using the tem peratures (or photon index ) prov ided 
bv lBrunschweiger et al l l|2009l ) and lTerada et all (|2008l ), re- 



spectively. We then extracted the hardness ratio from the 
spectra taking into account errors 3 (symmetric for all) by 
also reproducing the hottest and coldest spectra using the 
published errors for each source and computing the hard- 
ness. Before being able to add the BAT points to Fig. 
and [4] it is necessary to remove systematic differences be- 
tween the IBIS instrumental hardness and the BAT flux 
hardness. Fig. [S] shows the BAT hardness vs. the IBIS hard- 
ness for a sample of 13 IPs in common to both (we decided 
to excl ude IGR J00234— 614 1 due t o the extremely large er- 
rors in iBrunschweiger et al.l {2009;)). In red we display the 
one-to-one line where most of the data should sit in the ab- 
sence of systematic differences between the IBIS and BAT 
calibrations. However, it is easily noticeable that the BAT 
extrapolated hardness's are systematically harder than the 
IBIS ones. We compensate approximately for this by fitting 
a straight line through the datapoints and the origin (black 
line in Fig. [5]l . We then apply a correction to all the BAT 
points before plotting them in Fig. [2] [3] and [4] We note the 
large scatter about the fit, which might be caused by intrin- 
sic variability of the sources not being observed at the same 
time. This is also taken into account by folding the fit error 
into the BAT datapoints as well. 

All three plots show evident signs of correlations with 
hardness ratio when considering IPs alone. In particular, 
when considering IPs with Pspin/Porb < 0.1 (systems which 
are on their way to equilibrium at Pspin/Porb = 0.1 and are 
therefore exclusively in the propeller stage), then the corre- 
lations become even more evident. In order to obtain a sig- 
nificance for these correlations observed for the IPs thought 
to be in the propellor stage, we performed a Spearman rank 
test using a Monte-Carlo scheme. At first we decided to test 
the IBIS observations only, as we believe these are the mea- 
surements with lowest systematic errors. For example, in 



3 We note that normalisation constants are not required when 
inferring hardness ratios from single model spectra 



Hard X-ray mCVs 9 



— Data 

-- - Reference y-x 

— Fit y=0.58x 


















+ 



°0 0.5 1 1.5 

BAT Hardness 



error in the IBIS response, to any small spectral variabil- 
ity intrinsic to the observed systems. The addition of this 
intrinsic dispersion lowers the reduced chi-squared to unity 
and results in more conservative errors on the fit param- 
eters. In Fig(8] we display the contour plots for our linear 
fit in log-log space in the top panel, and the fit itself on 
the bottom. Note that the contours represent lines of a —1, 
2, 3, 4, 5. Again one can see that a simple constant value 
straight line fit is not consistent with the data. The result- 
ing equation to the fit can be expressed as 30-60/17-30 keV 
= ((Pspin/Porb) ~ 0.0259)-°- 21±0 05 - l0°- 09±0 - 03 , and may 
prove useful for modeling these systems and observations in 
the future. However at this stage, this is a purely empirical 
model. 



Figure 5. BAT extrapolated flux hardness vs. IBIS measured 
hardness. In blue is the data. In red we show a one to one line for 
reference. The cross-calibration fit is displayed in black. We note 
that both BAT and IBIS errors were taken into account when 
fitting. 



order to test if the correlation between hardness and P ap i„ 
is significant, we created 100,000 mock data sets containing 
the same number of points but shuffling the P ap in values 
randomly each time. Moreover, in order to take the hard- 
ness uncertainties into account, we replaced each hardness 
value with a random variable drawn from a normal distri- 
bution whose mean is equal to the observed hardness and 
whose standard deviation is equal to the error on the obser- 
vation. We then calculated the Spearman rank coefficients, 
p. The distributions of the coefficients for all P ap i n , Porb and 
Papin/ Port are shown in Fig. [6] Also displayed in each panel 
is the significance for the calculated Spearman rank value of 
the real data. These are 3.48cr, 2.74cr and 3.47cr for P 3p i„, 
Porb and Ps P in/Porb respectively. 

The fact that the correlation is apparent in all three 
plots is somewhat expected, since P sp in and P or b are not 
independen t, but expected to evolve together (jNorton et al.l 
120081 |2004| ) . At this stage we perform the same exercise as 
for Fig. H3 but this time including the additional IPs above 
the period gap observed by BAT only, and AE Aqr as ob- 
served by HXD. The results for this simulation are presented 
in FigEl 

As a final step we decided to extend our analysis further 
for the observed correlation in P sp in / Porb vs. hardness, given 
that, from an evolutionary p erspective, it is expected to be 
the most relevant parameter (|Norton et al.| [2008 . 2004j). We 
produce a linear fit in log-log space to the P sp in/ Porb vs. 
hardness plot in Fig. U for the hard X-ray selected mCVs 
used in this work 4 . When this is done, a non acceptable 
value of 4.5 is obtained for a reduced y 2 . H o wever , fol- 
lowing a simila r proc edure to iTremaine et al.l (|2002T ) and 
iMcHardv etHI l|2006l ). we introduce a small intrinsic dis- 
persion of 0.009 to all our datapoints in linear space. This 
corresponds to ~ 1% for the soft IPs in Fig. |3] and « 0.5% 
to the hardest IPs. There may be many reasons why such 
a small error might be introduced, ranging from a ~ 1% 



4 We have tried various polynomials but these all worsened the 
fit 



6 DISCUSSION 

IBIS has so far detected 32 CVs (23 spatially correlated with 
known CVs + 9 new, optically confirmed discoveries). The 
majority are intermediate polars, but IBIS has also detected 
the bright dwarf nova SS Cyg and a few polars. This sample 
is an ex tension of the previo usly presented sample of IBIS 
CVs by iBarlow etTal] l|2006h . which also showed that the 
spectral characteristics of these objects in the 20-100 keV 
range are actually quite similar and compare well with pre- 
vious high-energy spectral fits. Moreover BAT has observed 
a similar number of objects (with many in common), which 
shows that modern hard X-ray surveys are able to consis- 
tently observe mCVs at higher energies than before. 

The high incidence of IPs in our sample is not un- 
expected. Kuiipcrs & Pringlc (|l982h and many others, in- 
cluding IWarnerl ( 20031 ) , have suggested that the lower levels 
of hard X-rays/soft gamma-rays emission from polars may 
well be related to the low accretion rate and stronger mag- 
netic fields compared to IPs. It has also been suggested that 
the strong magnetic fields in polars a re able to produce a 
more "blobby" flow (|Frank et al.ll2002l ). These high density 
"blobs" are then able to penetrate deeper within the post- 
shock region and contribute more to the blackbody com- 
ponent of the broadband X-ray spectrum, and less to the 
bremsstrahlung component, making the broadband X-ray 
spectrum of IPs harder, and hence more luminous in the 
hard X-rays. 

Two out of four known APs are observed in our sam- 
ple. The remaining two happen to be i n a region of low IBIS 
exposure and are likely more distant. ISchwarz et alj (120051 ) 
have already shown that BY Cam, one of the observed APs, 
has different properties from those of normal polars, in par- 
ticular having a much higher accretion rate. We therefore 
tentatively conclude that IBIS has not yet seen the two 
missing APs due to their distance/exposure. Our sample 
only includes 2 definite synchronous polars, and we do not 
expect many of these systems to be observed in the future 
with higher sensitivities above 17 keV. 

Of the many observational characteristics presented 
here, one feature in the P sp in — P rb pl ane has stood out since 
the first study bv lBarlow" et all (2006) : a very low number of 
IPs below the period gap are detected with hard X-ray tele- 
scopes. The only exception in our study is the very nearby 
IP EX H ya. We can compare this r esult to the theoretical 
models of iNorton et all (|2008l. 120041 ). It is believed that the 
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Figure 6. Results from the Monte-Carlo simulation for estimating the correlation significances of the IBIS IPs only. The distributions 
display the calculated p coefficients for our mock datasets. We display with an arrow the calculated coefficient for the real set. 






Figure 7. Results from the Monte-Carlo simulation for estimating the correlation significances of the IBIS, BAT and HXD IPs. The 
distributions display the calculated p coefficients for our mock datasets. We display with an arrow the calculated coefficient for the real 
set. 



IPs below the period gap have accretion flows which are 
very different to the IPs above the period gap. In particular 
systems below the period gap with high P 3p i„ / ' P or b such as 
EX Hya display ring- like accretion. We would therefore not 
necessarily expect these systems to behave in the same way 
as the systems above the period gap, and, in particular, we 
would not necessarily expect them to emit such high energy 
ph otons. This is still a n open question, and, as mentioned 
by iNorton et al.l (|200&t ). only deeper hard X-ray exposures 
will reveal if this subclass of IPs displaying ring-like accre- 
tion is able to produce similar amounts of hard X-rays as 
those observed in disk-fed syste ms at longer orbital periods. 
It has already been mentioned l|Norton et al.|[2008h that as 
mCVs evolve through the period gap the magnetic field of 
the WD may be able to resurface when accretion stops. This 
can allow the system to synchronise, and we would then ex- 



pect a system jumping from P 3p i„ w 0.1P or i, above the pe- 
riod gap to Pspin ~ Porb below the period gap. Moreover 
we add to this that any system with P 3 pi n / P or b ^ 0.6 will 
never reach equilibr i um un til it reaches total synchronisation 
l|Norton et al.ll20Qg| , l2004h . 

All of the hard X-ray detected IP systems display 
Pspin < O.lPorb- This may be more observational evidence 
for different kinds of acc retion flows within the IP class, 
supporting the models of INorton et all (|2008l . |2004) . More 
evidence for these models comes from the non-detection of 
objects in any wavelength range within the synchronicity 
gap (a region above the period gap within P 3 pi n /P or b > 0.3 
and P 3 pin/Porb < 1). Such systems are predicted not to exist 
since IPs tend to evolve within the P 3 pi n — P rb plane towards 
their equilibrium spin rate at P 3p i„ « 0.1P or 6 above the pe- 
riod gap or at P 3p i n ~ 0.6P O rb below the gap. INorton et al.l 
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Figure 8. Top panel: Contours for a linear fit to the datapoints in the bottom panel. Lines display contours for a = 1, 2, 3, 4, 5. 
Bottom panel: IBIS IP hardness as a function of synchronicity. The equation resulting from the fit is 30-60/17-30 keV = {{Pspin/Porb) — 
0259) _0 21±0 - 05 — 10°' 09±0 - 03 



(2008, 2004) have predicted that low synchronisation mCVs 
have propeller accretion flows and are all trying to reach 
equilibrium moving towards P Bp in « 0.1P or b- This equilib- 
rium arises due to the WD trying to balance angular mo- 
mentum with the surrounding blobs. We believe this is the 
case for most hard X-ray selected IPs, as relatively few have 
yet been found above P sp in = 0.1P or b where the accretion 
flow is thought to take the form of a stream. 

Perhaps the most interesting result of this study is the 
discovery of a correlation between 30-60/17-30 keV hardness 
and spin/orbital parameters for IPs. No similar correlation 
has been reported before, probably because previously mea- 
sured X-ray hardness ratios of IPs were generally restricted 
to the range of approximately ~ 0.5 — 10 keV. Such ratios 
sample the lower end of the bremsstrahlung spectrum and 
the upper end of the blackbody spectrum, without fully mea- 
suring either component, while X-rays below ~ 1 — 2 keV are 
more easily affected by photoelectric absorption caused by 
the intervening accretion curtains. In contrast, we note that 
the spectral hardness variations we have measured in our 
hard X-ray detected IPs span the energy range ~ 17—60 keV 
and are only sampling the bremsstrahlung component of the 
spectrum. Therefore these observations tell us nothing about 
the relative contributions of the bremsstrahlung component 
emitted by the cooling plasma below the accretion shock 



and the blackbody component arising from the heated WD 
surface. Instead, they are directly sampling the relative con- 
tributions of the multi-temperature bremsstrahlung compo- 
nents that arise in the plasma below the shock front (given 
that the plasma cools as it settles towards the WD surface) . 

In order to explore a possible reason for the correlation 
between X-ray spectral hardness and spin-to-orbital period 
ratio in the hard X-ray detected IPs, we begin with the un- 
controversial idea that the WDs in IPs are mostly accretin g 
close to their equilibrium spin rates dNorton et al.l [2008). 
Hence, their spin-to-orbital period ratios are an indication 
of th eir magnetic field strength (see Fig. 6 of iNorton et al.l 
2008). Broadly speaking, for smaller P sp in/Porb, the surface 
magnetic field strength is smaller. So these systems will have 
smaller magnetospheric radii, and the material will attach 
onto field lines closer to the WD. This will give rise to a 
larger footprint area in those systems with smaller values of 
Pspin/Porb- That is to say, faster spinning WDs will have 
larger accretion footprints beneath a wide accretion cur- 
tain. Evidence for this also comes from the observed double- 
peaked pulse profiles observed in fast spinnin g WD (hence 
having a small magnetic field) as described in lNorton et al.l 
(1999), resulting from the optical depths across and along 
the accretion curtains as the WD rotates. Although the 
conventional theory of accretion shocks suggests that wider 
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shocks will be taller, we here investigate the consequences if 
the opposite is true. 

By spreading the material over a larger area beneath a 
shock of relatively low height, we suggest that the resulting 
bremsstrahlung X-ray emission may have a harder spectrum, 
possibly because the accretion is closer to the WD surface 
and there is less distance for the plasma to travel as it cools 
within the post-shock region towards the WD surface and so 
there is less contribution from cooler bremsstrahlung com- 
ponents. In contrast, the systems with a relatively slowly 
spinning WD have a larger Pspin/Porb value, so their mag- 
netic field strength is larger, their magnetospheric radius is 
larger, and their accretion footprint is smaller. We suggest 
that in such cases, the X-ray emitting region may sit beneath 
a tall but narrow accretion curtain, and then this geometry 
gives rise to a softer bremsstrahlung spectrum, possibly be- 
cause the accretion shock is further from the WD surface 
and so there is a greater distance for the plasma to travel 
within the post-shock region and cool as it falls towards the 
WD surface. This interpretation also helps explain the low 
detection number in the hard X-ray domain of EX Hya-like 
systems below the period gap with high P sp in/P rb which 
are thought to display ring-like accretion. In these cases the 
magnetospheric radius extends to a very large distance from 
the WD implying a very small footprint area. If this then 
means a very tall shock height (in line with the interpreta- 
tion above), then the plasma will have a long distance over 
which to cool as it travels towards the surface, and the spec- 
trum will be dominated by softer photons. 

Cont rary to this interpreta tion are the ideas pre- 
sented bv lSuleimanov et al.l l|2005l ) and lBrunschweiger et al.l 
(2009), where the spectral hardness in the hard X-ray range 
could be used in order to obtain the WD mass. However 
it is not yet clear if this idea could explain the hardness 
correlations presented here. In particular if the hard X-ray 
spectral hardness could be used to infer the mass of the 
WD primary, then we would also exepect the WD mass to 
correlate with Pspin, Port and P sp in/Porb, and this is not 
the case. Moreover the best evidence supporting this idea 
is presented in Figure 6 from Suleimanov et all (120051 ) and 
Figure 4 from iBrunschweiger et al. ( 20091 ). where the esti- 
mated WD masses inferred from the hard X-rays are plotted 
against the estimated WD masses inferred from other meth- 
ods. It is clear by looking at the plot that large scatter is 
present and that the error bars are quite large. We thus 
believe that no firm conclusion on the connection between 
WD masses and spectral hardness can be obtained at the 
moment. 

As far as we are aware, some det ailed modelling of ac - 
cret ion shocks ha s been performed bv lCanalle et al.1 (|2005l ) 
and lSaxton et all (120071 ). however their model X-ray spectra 
obtained by the simulations only extends to 10 keV, which is 
not high enough for the data considered here. We point out 
that no one has yet modelled whether the multi-temperature 
bremsstrahlung spectrum is different for a wide, low accre- 
tion curtain compared with a tall, narrow accretion curtain, 
but we suggest this would be a worthwhile test to carry out. 

We recognise that our suggestion of wide/low accretion 
regions in the case of rapidly spinning white dwarfs, and 
narrow/tall accretion regions in the case of slowly spinning 
white dwarfs is contrary to the prediction of standard shock 
theory that wider accretion regions are taller. However, there 



is much in the detailed accretion flows of magnetic CVs that 
remains poorly understood and we put our suggestion for- 
ward in the spirit of advancing understanding of the subject. 
It then provides a possible way of understanding the hard- 
ness - period correlation observed here which would other- 
wise be unexplained. 



7 CONCLUSIONS 

This paper has presented a catalogue and analysis of a sam- 
ple of CVs detected in the hard X-ray range (> 17 keV) with 
IBIS, BAT and HXD. As with previously compiled high- 
energy samples of CVs, it is shown that most systems are 
magnetic. Moreover, some of the detected systems are very 
rare types of objects (2 APs). The sample is dominated by 
intermediate polars, with only 2 synchronous polars. This 
suggests the broadband X-ray/gamma-ray spectrum of IPs 
is harder and more luminous than that of polars. This could 
be the effect of accretion rate and magnetic field strength, 
where IPs have higher accretion and weaker magnetic fields 
relative to polars, depositing a somewhat smoother accre- 
tion stream onto the poles. By contrast to the polars, the 
accretion flow in IPs may therefore not bury itself deep 
within the post-shock region, producing a harder broadband 
X-ray /gamma-ray spectrum. 

We have shown that only IPs with P sp in/P rb < 0.1 
are consistently found by hard X-ray surveys. Moreover, 
we have examined the observational properties of mCVs 
in the P or b — Pspin plane. We find that the observations 
are consiste nt with the theoretical models of iNorton et al.l 
l|200Sl . |2004T ) for mCV evolution, where IPs tend to cluster 
at about P sp in/Porb ~ 0.1, and none have yet been observed 
in the hard X-ray regime above Pspin/Porb ~ 0.1. Also ob- 
served and predicted is the observation that a very low num- 
ber of IPs are found in any wavelength range within the 
synchronicity gap: a region between Pspin/Porb ~ 0.3 and 

Ps P in / Porb — 1- 

The paper has also presented the first observed cor- 
relations between 30-60 keV/17-30 keV hardness and Pspin, 
Porb and P ap i n I 'Porb- The correlations have been statistically 
tested using Monte Carlo simulations. 

In an attempt to explain our result we have suggested 
that hard X-ray selected IPs are spinning towards their equi- 
librium, so that their spin period is an indicator of magnetic 
field strength. This in turn will give rise to a short, but wide, 
post-shock region for fast spinning WDs (and therefore pos- 
sessing a relatively weak magnetic field) making their hard 
X-ray spectra harder. In contrast slowly spinning WDs will 
have a tall but narrow post-shock region (possessing a rela- 
tively high magnetic field), yielding a cooler bremsstrahlung 
component in the hard X-rays. 

All of the observations presented in this paper are con- 
sistent with mCV evolution models. It is very likely that 
hard X-ray missions will continue to increase this sample of 
mCVs, and it is also expected that more unidentified hard 
X-ray sources will be identified as IPs with more optical 
follow-ups. In particular, more observations will allow us to 
establish if any IPs are detected below the period gap and 
if any IPs will ever get detected in the IBIS energy range 
above Pspin /Porb ~ 0.1 in order to establish if the Norton 
accretion models are a plausible explanation to the observed 
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systems. This also implies that hard X-ray selected samples 
could have their own biases, however more analysis will have 
consequences on evolution studies of these exotic magnetic 
systems. 
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